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Abstract 

A theory is presented for the membrane junction separation induced by the adhesion between 
two biomimetic membranes that contain two different types of anchored junctions (receptor/ligand 
complexes). The analysis shows that several mechanisms contribute to the membrane junction 
separation. These mechanisms include (i) the height difference between type-1 and type-2 junctions 
is the main factor which drives the junction separation, (ii) when type-1 and type-2 junctions have 
different rigidities against stretch and compression, the "softer" junctions are the "favored" species, 
and the aggregation of the softer junction can occur, (iii) the elasticity of the membranes mediates 
a non-local interaction between the junctions, (iv) the thermally activated shape fluctuations of the 
membranes also contribute to the junction separation by inducing another non-local interaction 
between the junctions and renormalizing the binding energy of the junctions. The combined effect 
of these mechanisms is that when junction separation occurs, the system separates into two domains 
with different relative and total junction densities. 

PACS numbers: pacs numbers: 87.16.Dg, 68.05.-n, 64.60.-i 
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I. INTRODUCTION 



Adhesion of membranes is responsible for cell adhesion which plays an important role 
in embryological development, immune response, and the pathology of tumors 0. In 
many cases, membrane adhesion in biological systems is mediated by the specific attrac- 
tive interactions between complementary pairs of ligands and receptors which are an- 
chored in the membranes. At the same time, the adhesion between multi-component 
bio-membranes or biomimetic membranes is also intimately related to domain forma- 
tion. [Rl f|, |5], [7|, ^ |S], [IT], [nj When the membrane adhesion is mediated by the specific lock- 
and-key type of bonds between the anchored ligands and receptors, i.e., junctions, adhesion- 
induced lateral phase separations have been observed in many experiments in biomimetic 
systems |§, §], ||. Theoretical models and Monte Carlo simulations f| 0, || |9|, [[(J, [H] have 
also shown similar phase separation behavior in various systems. 

So far, studies on adhesion-induced lateral phase separation have focused on the case when 
the system has a single type of junctions. The presence of the glycol proteins anchored in the 
membranes (i.e., repellers), and the interplay between generic interactions (for example, Van 
der Waals, or electrostatic interactions) and specific ligand/receptor interactions are believed 
to enhance this phase separation. However, in biological systems membrane adhesion are 
often mediated by more than one type of junctions, and the adhesion-induced junction 
separation are believed to play an important role in some biological processes. For example, 
a key event governing a mature immune response when T lymphocytes interact with antigen- 
present cells is the formation of immunological synapses. An immunological synapse is a 
patch of membrane adhesion region between a T cell and an antigen-present cell, where the 
TCR/MHC-peptide complexes aggregate in the center with a LFA-l/ICAM-1 complexes rich 
region surrounds it. [[TJ, |H| Q Since a complete understanding of the physical mechanism 
behind this type of adhesion-induced multi-species membrane junction separation is still 
unavailable, in the present work I develop a theoretical model to study the equilibrium 
properties of such systems. 

This article is organized as follows. In section O, I discuss a coarse grained model for 
the adhesion of two membranes due to the formation of two types of junctions. To concen- 
trate on the effect of the differences between type-1 and type-2 junctions, the glycocalyx 
and the generic interactions between the membranes are not considered in this model. Fur- 
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thermore, I assume that the membranes are bound to each other due to the formation of 
the membrane junctions. Hence I will not discuss another interesting problem of the un- 
binding transition. An approximate solution of this model which neglects the fluctuations 
of membrane-membrane distance (the "hard membrane" solution) is studied in section [TQ|. 
This simplified solution already reveals several mechanisms which are important to the phase 
behavior of the system. For example, when type-1 and type-2 junctions have the same rigid- 
ity, membrane adhesion can induce a junction separation which is driven by the height 
difference of the junctions. In this situation the membranes separate into a type-1- junction- 
rich domain and a type-2-junction-rich domain. On the other hand, when the junctions have 
different rigidities and the height difference is not very large, membrane adhesion can induce 
an aggregation of the "softer" junctions, i.e., the membranes can separate into two domains 
which are both rich in the softer junctions. Thus, in general situations both mechanisms 
contribute to the adhesion-induced junction separation. When phase separation occurs, the 
system separates into two domains with different membrane-membrane distance because of 
the height mismatch of the junctions, and the total number of the softer junctions in the 
system is greater than the total number of the stiffer junctions. 

The fact that the hard membrane solution assumes that the membrane-membrane dis- 
tance is a constant has its drawback, too. An apparent artifact of the hard membrane 
solution is that when junction separation occurs, the total junction density, i.e., 0i + 02 
(4> a is the density of type-a junctions), is the same in domains which have different values 
of (0i — 02) /(0i + 02) (the relative densities of the junctions)! This artifact also shows 
that, in order to acquire a complete physical picture of the junction separation, the effect of 
non-constant membrane-membrane distance, and the thermally activated fluctuations of the 



membranes and junction densities should be taken into account. Therefore in section [TV] I 
study the effects of these fluctuations. The fluctuation analysis shows that, first, the ther- 
mally activated membrane fluctuations renormalize the chemical potentials of the junctions, 
and effectively reduce the binding energies of the junctions. This chemical potential renor- 
malization is less significant for the softer junctions because they allow the membranes more 
freedom to move. Second, the fluctuation analysis also reveals nonlocal interactions between 
the junctions which are mediated by the membrane elasticity and thermally activated fluc- 
tuations of the junction densities and membrane-membrane distance . These interactions 
are not included in the simple physical picture provided by the hard membrane solution. 
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As a result of these effects, when junction separation occurs, domains with different values 
of (4>i — 4>2)/{<j)\ + <t>2) also have different values of 4>i + 4>2- The fluctuation analysis also 
shows that, when the hard membrane solution of the junction densities are small, or when 
the junctions are very short or very soft, the membrane fluctuations are sufficiently large 
such that the present analysis cannot provide the complete physical picture for the system. 
This criterion shows under what conditions one needs a numerical simulation of the model 
to provide a better picture of the physics in this system. Section |V] summarizes this work. 
The Appendix discusses the details of the fluctuation analysis around the hard membrane 
solution. 



II. THE MODEL 

To focus on the physics of adhesion-induced phase separation, I will not discuss the 
binding/ unbinding transition but only consider the case when the membranes are bound to 
each other due to the presence of the junctions. The system is shown schematically in Fig. 1. 
The heights of the membranes measured from the reference plane (i.e., the x?/-plane) are 
denoted as Zi(r) and 22 (r), respectively, where r = (x, y) is a two-dimensional planar vector. 
There are two types of anchored receptors in membrane 1, and two types of anchored ligands 
in membrane 2. Type-a receptors (a is 1 or 2) form specific lock-and-key complexes with 
type-a ligands, these are the junctions which mediate the membrane adhesion. The density 
of type-a junctions at r is 4> a (r), and the densities of free type-a receptors and ligands at r 
are denoted by ipRai?) and i>La{ v )i respectively. The binding energy of a type-a junction is 
denoted by Ebu- 

The effective Hamiltonian of the system can be written as 

H= \ Sr \\ ( y2/ ^) 2 + 1 ( v/ ^ 2 + ^ Y Mr) {k{r) - ha)2 - £ ■ (1) 

The energy unit is chosen to be fcgT. Here h(r) = -Zi(r) — -22(1") is the membrane-membrane 
distance at r. The first and second terms on the right hand side are the bending elastic energy 
and surface tension of the membranes, k is related to the bending moduli of the membranes 



by k = k\K2/{k\ + K2) jnj , and 7 is related to the surface tension of the membranes by 



7 = 7172/ (71 + 72) In this simple model it is assumed that k and 7 are independent of 



the densities of the receptors and ligands anchored in the membranes. I also assume that in 
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the presence of a type-a junction, the interaction energy between the membranes acquires 
a minimum at h — h a (the natural height of a type-a junction), and the coupling term 
J2a=i ^0o( r ) (M r ) — ^«) 2 comes from the Taylor expansion around this minimum. Here 
X a is the rigidity of a type-a junction against stretch or compression. The last term on 
the right hand side is the binding energy between the receptors and the ligands. To focus 
on the effect of adhesion-induced interactions, I have neglected all the direct interactions 
between the junctions, receptors, and ligands. The nonspecific interactions between the 
membranes are also neglected. For simplicity, from now on I further choose the unit length 
in the xy-plane to be y/a, where a is the in-plane size of an inclusion, and the unit length 
in the z-direction is chosen to be \Jcl/k = Iq- Thus the Hamiltonian of the system can be 
expressed in the non-dimensional form, 

H = Jd\l l - (V 2 Mr)) 2 + V - (VMr)) 2 + £ ^0 a (r) (h(r) - h a f - £ <P a E Ba ) , (2) 

where Y = 7Z 2 , is the dimensionless surface tension, A a = \ a ll is the dimensionless junction 
rigidity, and all in-plane lengths and heights are scaled by yfa and \fa/~K = l , respectively. 

The effective interaction free energy between the junctions due to the membrane-junction 
coupling is obtained by integrating over h(r), 

F c [<f> a ] = -\n(jD[h]e- H ^y (3) 

Thus in the spirit of density functional theory, the total free energy of the system is provided 
by 

F = F C + F„ (4) 

where 

2 r 

F s = J2 rfV (0 Q (r)(ln0 a - 1) +^ Ra (r)(\mlj Ra - 1) + ^(r) (In ^ - 1)) (5) 

a=l 

is the contribution from the entropy of the junctions, receptors and ligands. Here I have 
assumed that ^ tt <C 1, ipR a <S 1, and i^La ^ 1- In principle, once F c is calculated, the 
equilibrium distribution of the junction density is determined by minimizing the total free 
energy of the system under the constraint that the total numbers of the receptors and ligands 
in the system are fixed, i.e, 

J d 2 r{<f) a {r) +4> Ra (r)} = N Ra , 

J rf 2 r{0 a (r) + ^ La (r)} = N La , (6) 
5 



here Nji a and Nl q are the total number of type-a receptors and ligands in each membranes 
when the membranes are completely detached. 



III. "HARD MEMBRANE" SOLUTION 

Since the integral in Eq. (^) cannot be carried out exactly, in this section I discuss an 
approximate solution in which cf) a and h(r) are independent of r. In this approximation, F c 
can be easily calculated by looking for the saddle point in the integrand. This is equivalent 
to neglecting the fluctuations of the membrane-membrane distance, therefore I call this 
mean-field approximate solution the "hard membrane" solution. To simplify the notation, 
I define A± = Aj ± A 2 , <f>± = (0i ± <f>2)/2, and let hi = ho — Ah, hi = h + A^. Thus the 
hard membrane solution of the membrane-membrane distance can be expressed by 

A_0 + + A + 0_ 
" = /l °- A + + + A_0_ Afe 

= h - l M . (7) 

Notice that Im depends on the junction densities. After substituting h back to the Hamil- 
tonian, the effective interaction free energy between the junctions, F c , can be expressed by 
its saddle-point value 

d 2 r j A +^+ + A -^- (4 + A 2 h ) + (A_0 + + A + 0_) l M A h - + j , (8) 

where E± = Ebi ± Eb2- It is clear that there is an interaction between the junctions due to 
the membrane adhesion. To minimize the total free energy under the constraints in Eq. (^j), 
it is convenient to work in the grand canonical ensemble and define the free energy G of the 
system under constant chemical potentials, 

G = F c + F s -Ys^Ra d 2 r((j) a + 4j Ra ) -Y^^La d 2 r(4> a + tp La ). (9) 

The chemical potentials, /xr q , fiia are determined by fixing the total number of receptors 
and ligands in the system. However, for convenience I will proceed the discussion in the 
grand canonical ensemble. After some straightforward algebra, G is expressed as 

G = Jd 2 r 0+ {g(4>) + 2 (In 0+ - 1) - //+} + G^ (10) 
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where 

g(4>) = -^A + ^^ + (l + 0)ln(l + 0) + (l-0)ln(l-0)- / i_0 

= /(0)-^_0, (11) 

A 2 

A = A_/A+, 4> = 4>-/4>+, and fj,± = (/i R1 + fi L1 ) ± (fi R2 + /x L2 ) + #bi ± E B2 - ^A ± . 

includes terms which only depend on ip Ra and -0La, they are decoupled from the other 
terms, hence from now on I neglect Gf. From Eq. fllOD, it is clear that in the hard membrane 
solution the phase behavior of the junctions is governed by A^A + , A, and n±. Minimizing 
g{4>) leads to the equilibrium value of 0, and later I will show that there can be a phase 
separation in 0. On the other hand, from Eq. (|10l), 0+ satisfies 

0+ = exp - -g((f>)) . (12) 



X T 

Because in equilibrium <p is determined by minimizing g(<f>), g{4>) takes single value even 
when there is a phase separation in 0. Therefore in the hard membrane solution + is 
single-valued even when the system separates into two domains with different values of 
0! In the next section I will show that, when the effects of fluctuations around the hard 
membrane solution are taken into account, the analysis reveals a renormalization of the 
binding energy of the junctions and (nonlocal) interactions between the junctions which 
are mediated by the membrane elasticities and thermally activated membrane fluctuations. 
As a result, the true equilibrium solution of + is not single- valued in the regime where 
junction separation happens. Thus, the fact that 0+ is single-valued in the hard membrane 
solution is an artifact of the approximation which assumes constant junction densities and 
membrane-membrane distance. 

Now I discuss the hard membrane solution of 0. To emphasize different roles played by 
A|A + and A, I begin the discussion with the special case when A = 0, i.e., when both types 
of junctions have the same rigidities. In this case the important parameter of the theory is 
A|A + , and g(4>) has a very simple form 

9(4) = + (1 + 4>) In (1 + 0) + (1 - 0) In (1 - 0) - //_0. (13) 



This form is exactly the same as the Flory-Huggins theory for binary mixtures |L7j, where 
phase separation occurs when A|A + > 2 and the phase coexistence curve is a straight line 
at fi_ = 0. This phase coexistence curve ends at a critical point = 0, A^A + = 2. The 



physics in this special case A = is clear: the difference in junction height drives a junction 
separation, and this separation only occurs when the factor A|A + , a combination of junction 
height difference and junction rigidity, is sufficiently large. On the phase coexistence curve, 
the system separates into 0!-rich and 2 -rich domains, and the system is symmetric under 
0^-0. 

Next I discuss the more general case A 7^ 0, i.e., the junctions have different rigidities. 
Fig. 2 shows the shape of g(4>) with different values of /i_ when A = 0.2 and Aj;A + = 1.998. 
Notice that this is the case when A|A + < 2, i.e., there is no junction separation if A = 0. 
Nevertheless, Fig. 2 clearly shows that g(4>) has two local minimum, both at negative 0, and 
phase coexistence occurs when /i_ pa —0.4045. Since this is the case when A > 0, i.e., type-2 
junctions are "softer" than type-1 junctions, double minimum at = 0i — 02 < means 
that the softer junctions tend to aggregate, when phase coexistence occurs a domain with 
high 02 — 0i coexists with a domain with small 02 — 0i. The density of type-2 junctions is 
higher than the density of type-1 junctions in both domains. 

Another effect of nonzero A can be seen in Fig. 3 , where g{4>) for different values of 
A is shown at A^A + = 2.04 > 2. It shows that g(4>) is symmetric in when A = but 
asymmetric in for nonzero A, i.e., the symmetry under — > —0 no longer exists when the 
junctions have different rigidities. Comparing to A = case, in the case when A > 0, the 
minima of g(0) are shifted towards smaller values, i.e., the softer junctions are easier to be 
formed. Notice that different from the example in Fig. 2, in Fig. 3 when phase coexistence 
occurs the membranes separate into 0x-rich and 2 -rich domains, but the softer junctions 
(in this case type-2 junctions) are the "favored" species, i.e., the total number of the softer 
junctions in the system is greater than the total number of the stiffer junctions. From these 
two examples of nonzero A, I conclude that in general the experimentally observed junction 
separation induced by membrane adhesion is actually a result of the combined effect of the 
aggregation of softer junctions and the separation of the junctions due to the mismatch of 
junction heights. 

In the neighborhood of A^A + = 2, A = 0, the equilibrium value of is small compared to 
unity, therefore the phase diagram of the system in this regime can be studied by expanding 
g(0) around = 0, 

g(</>) = r 2 2 + r 3 3 + r 4 4 - /i_0 + const. + C(0 5 ), (14) 



8 



r 2 = l-fA + (l-Af , 

rs = ^A + A(l-A 2 ) 2 , 

, 4 = I-fA + A*(l-Af, 

= + ^A + A (2 - A 2 ) , 

A 2 

const. = — ^A+A 2 , (15) 

and 0(<f) 5 ) is the contribution from terms of order 5 and higher. The phase diagram in 
the neighborhood of A|A+ = 2, A = is plotted schematically in Fig. 4. where the phase 
coexistence curve for A = ends at a critical point A^A + = 2, //_ = 0, and the end 
points of the phase coexistence curves for A ^ occurs at the triple root of dg/d<p = 0. 
Straightforward calculation leads to the position of the end points of the phase coexistence 
curves at 

A 2 A+ = 2(l-9A 2 /4) + 0(A 4 ), 

/i_ = -2A + C(A 3 ). (16) 

This shows how the smallest value of A 2 A + above which phase separation can occur decreases 
as the difference of junction rigidities increases. The phase coexistence curves move towards 

A 2 A 2 

the A = phase boundary as the value of -^^+ increases. This is because as -^-A + increases, 
the effect of junction height mismatch becomes more important, and the difference in the 
junction rigidities becomes less important. 

Although the hard membrane solution is a very simplified analysis of the model, it nev- 
ertheless reveals interesting physics of the junction separation due to membrane adhesion. 
First of all, the height difference between different types of junctions is not the only factor 
which is important for the junction distribution. It is only when the rigidities of both types 
of junctions are the same that the difference in the junction height is the most important 
factor in the junction separation. If the system consists of junctions with different rigidities, 
the softer junctions are more favored, and the junction separation is a result of the interplay 
between the aggregation of the softer junctions and the separation of the junctions due to 
the height difference. Therefore the smallest value of A|A + above which phase separation 
can occur is smaller for systems with larger |A|. However, the approximations in the hard 



9 



membrane solution lead to the surprising result that 0+ is the same in both domains when 
junction separation occurs. This indicates that analysis which includes the elasticity of the 
membranes and the thermal fluctuations around the hard membrane solution is very much 
needed for a full understanding of the nature of the membrane adhesion-induced interactions 
between the junctions. 

IV. BEYOND "HARD MEMBRANE" SOLUTION 

As mentioned in the previous section, the hard membrane solution neglects the effects of 
nonuniform membrane-membrane distance and junction densities, therefore fails to take the 
effects of membrane-mediated nonlocal interactions between the junctions into account. The 
result is reflected in the unrealistic solution of single valued + in both domains when phase 
coexistence occurs. To study these membrane-mediated effects, in this section I include the 
fluctuations of membrane- membrane distance and junction densities by expanding the free 
energy of the system around the hard membrane solution. In the following I denote the true 
membrane-membrane distance as 



Here 51, 5<j) a are the deviations of the true values of h and 4> a from their hard membrane 
solutions, (p a M and hu are the hard membrane solution of Q (r) and h(r), respectively. In 
this expansion, the coarse grained Hamiltonian can be expressed as 



h(r) = h + l M + 51(t) = h M + Sl{r), 



(17) 



and the densities of the junctions are expressed by 



0a = (paM + <50(r). 



(18) 



H = H M + H + Hi + H^, 



(19) 



where H M is H(h M , <Pim, 4>2m), 




includes terms which are bilinear in 81 and 5<j), 




(21) 
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is the nonlinear coupling between 51 and 5(j) a , and includes terms which are linear in 

5(j) a . 

First I discuss the contribution from H , i.e., the Gaussian fluctuations around the hard 
membrane solution. In this Gaussian approximation, F c has acquired two correction terms 
which can be expressed by 

1 V In 2?F T \l M 5m + (q) + A h (5m_(q)| 2 

2, ? 4 + r 9 2 + m + Y q i + Tq 2 + m + ' 1 ' 

for convenience I have defined m± = K\(j)iM ± A 2 2 a/, an d 5m± = Ai5<pi ± A 2 <50 2 - The 
first term is independent of 5<fr a , therefore I neglect it in the rest of the discussion. The 
second term is a membrane- mediated nonlocal interaction between the junctions. This 
interaction has two characteristic lengths: m + ^ 4 is the distance it takes for a perturbation 
in membrane-membrane distance to relax back to its hard membrane solution due to the 
membrane bending rigidity, another length is r _1//2 , for lengths greater than T~ 1//2 the 
elasticity of the membrane is dominated by the surface tension of the membrane, and the 
contribution from the bending rigidity is negligible. In the rest of this article, I focus on 
the case when T < ^/m + , in which the membrane bending rigidity is the dominant effect 
which drives a perturbation in h back to thus the contribution from surface tension 
of the membranes is negligible. To understand the nature of the nonlocal interaction, it is 
convenient to transform the second term to real space. Calculations in the Appendix show 
that, when T < ^m + , the second term in the real space has the form which is derived in 
Eq. (H) 

1 _ iM A^r) - (l + ^) A 2 50 2 (r) 



/ { f v I d 2 r' A ^ G(|r-r>y 4 ) 



X 



— W<Mr')- 1 + — A 2 50 2 (r') 



(23) 



where G{x) is a MeijerG function fl8| . G(x) is vanishingly small for x > 5. Eq. (|23|) shows 
that this membrane-mediated interaction is attractive between junctions of the same type, 
and repulsive between junctions of different types. This interaction is short-ranged with a 
characteristic length m + 1,/4 . Also notice that the contribution from H\ is proportional to 
A^, i.e., there is no membrane-mediated interactions in the level of Gaussian approximations 
when the junctions have the same height. The physical picture of this interaction can be seen 
from Fig. 5. A small perturbation of the junction density from the hard membrane solution 
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induces a deviation of membrane-membrane distance from Hm, and there is a membrane 
bending energy associated with any given distribution of non-uniform membrane-membrane 
distance. To reduce the bending energy, a region with positive 5(f)^ 2 ) attracts a region with 
positive 50i(2) i n order to reduce the elastic energy cost of a "pit" or a "bump" between 
these two regions due to the non-uniform h(v). Similarly, a region with positive S(f>i repels a 
region with positive 5(j) 2 , in order to reduce the bending energy cost due to the high curvature 
configuration between these two regions. This also explains the fact that these interactions 
vanish when both types of junctions have the same height, i.e., = 0. A similar kind of 
membrane-mediated nonlocal interaction between the junctions is discussed in the celebrated 



article by Bruinsma, Goulian, and Pincus [19], where in their "van der Waals regime", the 



competition between the potential minimum due to the van der Waals interaction between 
the membranes and another potential minimum due to the stiff membrane junctions results 
in a strong interaction between the junctions. Although there is only one type of junctions in 
the system discussed in Ref. [|l9j, the interaction between the junctions in Ref. [jnj and the 



present case share the same physical mechanism, i.e., the bending elasticity of the membranes 
mediates this interaction. 

Another type of nonlocal interactions between the junctions can be studied by consid- 
ering the effect of nonlinear couplings between 51 and 5cj) a . This is done by including the 
contributions form Hi perturbatively to one loop order. The resulting effective interaction 
free energy between the junctions, F c , now has the form 

F c = F M + Fq + F loop + (24) 

where Fm is the hard membrane solution of F c , Fq is the contribution from terms which are 
bilinear in 51 and 5(f) a , and Fi oop is the contribution from the nonlinear couplings between 
51 and 5(p a to one loop order. The details of the calculations for F\ oov is discussed in the 
Appendix. When T < v /m+, the result (up to terms quadratic in 5(p a ) is provided by 
Eq. ©, Eq. (g7|), and Eq. ©, 

Fi oop = - 1 / rf 2 r(A 1 50 1 + K 2 5<p 2 ) 

1 f d 2 q 1 



lb^/m + J (2it) z g 4 + 4m + 
Here the first term is a "renormalization" of the chemical potentials of the junctions due to 
membrane fluctuations. This term effectively reduces the binding energies of the junctions. 
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The fact that the renormalization of the chemical potential for the softer junctions is less 
significant compared to that for the stiffer junctions is because the membrane fluctuations 
are energetically less costly for the softer junctions. The second term is a fluctuation-induced 
nonlocal interaction between the junctions, and higher order terms are neglected. Notice 
that, as discussed in the Appendix, the second term in Eq. is actually an approximate 
form of the much more complicated true result, it provides the correct asymptotic behavior 
of the true result at large and small q limits in the case when V < ^m + . Similar to 
the case of Gaussian approximation, when the fluctuation-induced interaction between the 
junctions is expressed in real space, one finds that the interaction between the junctions is 
nonlocal, short-ranged, and has a characteristic length on the order of m^ 1 ^ 4 . Since Fi oop is 
non- vanishing even when A^ = 0, it is clear that the thermal fluctuations of the membrane- 
membrane distance is the mechanism which induces the nonlocal interactions between the 
junctions in Fi oop . This is similar but not the same as the interaction between the junctions 



in the "Helfrich regime" discussed in Ref. |T{|. In Ref. []T9|| , the interaction between the 
junctions in the Helfrich regime comes from the collisions between the membranes. Here in 
the one-loop calculation the interaction between the junctions comes from the fluctuations 
of the membrane-membrane distance around the hard membrane solution, the effect of 
membrane collisions is not included. 

When the fluctuations around the hard membrane solution are taken into account to 
one-loop order, the total free energy of the system to second order in S(j) a can be expressed 
by 

d 2 q 



J2-—\suq)\ : 

a=1 ^<PaM 



(2vr) 2 

(2tt)V + ™+ 
' ; " <H] ' |A 1 ^ 1 ( g )+A 2 50 2 (g)| 2 



-— )A 1 6cf> 1 (<l)-(l + — WW 
m_i_ / V m+ / 



IQ^fm^J (27r) 2 g 4 + 4m^ 
+ - 1 / d 2 r(A 1 5(f) 1 + A 2 5<p 2 ). (26) 



Here the first term on the right hand side is the hard membrane solution, the second term 
comes from the entropy of the junctions, the third term is the nonlocal interaction between 
the junctions due to Gaussian fluctuations, the fourth and the fifth terms come from the 
nonlinear couplings between SI and 5cj) a . Notice that the contribution from does not 
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appear in the total free energy of the system, it cancels with the linear terms in the expansion 
of the entropy of the junctions. This is because <p a M minimizes the hard-membrane free 
energy, therefore in the expansion around the hard-membrane solution terms which are 
linear in 50 cancel with each other. The contribution from one-loop calculation, however, 
includes terms which are linear in 50 a because they come from the nonlinear couplings 
between 50 Q and 51. An important consequence of the presence of these terms is that in 
general the equilibrium values of 50i and 502 are nonzero due to the membrane fluctuations. 
Therefore when a phase separation occurs, the values of + 02 are different in domains 
with different values of 0. 

To discuss the correction of <fi a and h(r) due to the thermally activated fluctuations, it is 
convenient to express Eq. (p6|) as 



72 

F = F M + J d 2 r (<5M0i(r) + W^(r)) + £ / —^M a/3 (g) 50^)50^), ( 2 ?) 

af3 v / 

where 



A A 
6 H a 



+ 



1 A? / m\ 2 a2 1 A'i 



hM q + m + \ m + J lQ^/m + g 4 + Am. 

1 A2 / m\ 2 a2 1 A 2 
M 22 (q) = — 1 + A 2 2 - — 

20 2A f q* + m+\ m+J 2 16 v /^g 4 + 4m H 



Now S(j) a (q) can be expressed by 8fj, a and M a/ g, 

50x (g) + 50 2 (g) = 5(g) x 2d ~ 1 M {(M 22 - M 2 i) 6 m + (M n - M 12 ) 5/i 2 } , 

50i(g) - 50 2 (g) = 5(g) x {(M 22 + M 21 ) 5^ - (M n + M 12 ) 5/i 2 } , (29) 

where det M = M11M22 — M12M21. This rather complicated expression shows that, besides 
//_, A, and Aj l A + , the answer to the question of which domain has higher total junction 
density when the phase coexistence occurs also depends on the values of A 2 and a M (to 
determine <fr a M, one needs to know the value of //+)! In this article I shall not discuss the 
details of the values of 0i + 02 for different given parameters in the theory, but simply 
comment that when det M is positive, the phase diagram of the hard membrane solution is 
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not modified by the thermal fluctuations. However, when det M < 0, the hard membrane 
solution is not stable at any finite temperature. The result in Eq. ( p9|) also provides some 
criteria for the current analysis. For example, when the fluctuations are large, the deviation 
from hard membrane solution can no longer be treated by perturbation theory. This is 
true when S(f) a /(p a ~ 0{\). Since 5(f> a becomes large for small detM, which occurs at 
small m + = Ai0ia/ + A202M, I conclude that the perturbation theory breaks down at 
small junction densities. Finally, I point out that the collisions between the membranes are 
also neglected in the present analysis, this approximation is valid when the fluctuations of 
membrane-membrane distance is not large, i.e., when 



<(Siy> ( r d 2 q 1 \ 1/2 1 1 



/ q 1 

h M \l (27r)V + m J ^ - ( } ' { ) 

In the regime where rnj^hu = (M4>im + A 2 02m) 1 ^ 4 ^a/ < i- e -; when the junction 

densities are small, or the when junctions are very soft, or when the junctions are very 
"short", the contributions from membrane collisions should be taken into account for a 
complete analysis of this system. Thus, when the membrane fluctuations or the membrane 
collisions become important, numerical simulations |20|| or other methods which take the full 
membrane fluctuations into account should be applied to study the physics of this system. 



V. SUMMARY 



I have discussed the phase separation of multiple species membrane junctions induced 
by membrane-membrane adhesion with a continuum theory. In the hard membrane ap- 
proximation where the membrane-membrane distance and junction densities are assumed 
to be constants, we find that A|A + and A are the important parameters which governs the 
junction separation. When A = 0, both types of junctions have the same rigidity, and the 
junction separation is driven by the height difference of the junctions. Under this condition 
the junction separation is very similar to the Flory-Huggins theory for a binary mixture. 
Phase separation occurs when A|A + > 2 and |U_ = 0, the phase coexistence curve ends at 
a critical point fi- = 0, A|A + = 2. When A ^ 0, the junctions have different rigidities, 
and the softer junctions are easier to form than the stiffer junctions. Therefore the softer 
junctions have a tendency to aggregate. In this more general case, the height difference and 
the junction rigidities difference both drive the phase separation, thus the phase separation 
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can occur at A|A + < 2. 

The Gaussian fluctuations around the hard membrane solution reveals a membrane- 
mediated nonlocal interaction between the junctions. This interaction is short-ranged, which 
decays with a characteristic length (A + + a/ + A_</>_m) -1 / 4 , it is attractive between the same 
type of junctions, but repulsive between different types of junctions. The strength of this in- 
teraction is proportional to A|, and it is due to the membrane bending energy cost between 
regions with different junction densities. Perturbation theory to one-loop order shows other 
effects of thermal fluctuations, this includes a renormalization of the chemical potential of 
the junctions which effectively reduces the binding energies of the junctions, and a nonlocal 
interaction between the junctions which is independent of A/j. The fact that the contribu- 
tion from one- loop calculation is non- vanishing even when junctions of type-1 and type-2 
have the same height indicates that this contribution is a result of thermal fluctuations of 
the membranes. Hence it is non-vanishing at all finite temperatures. When the contribution 
from one-loop calculation becomes very large, the hard membrane solution is qualitatively 
incorrect, and the effects of thermal fluctuations is a dominant factor. This can occur at 
very low junction densities. The Gaussian fluctuations of the membrane-membrane distance 
also provide another limit of the present analysis: the mean squared fluctuations of the 
membrane-membrane distance should be small compare to Hm- As a result, the analysis in 
this article does not provide the complete physical picture of the system for very soft or very 
short junctions, either. 

In summary, mean field and fluctuation analysis of a simple coarse grained model for 
adhesion-induced phase separation of multiple species of membrane junctions is studied 
in this article. This model shows rich behaviors which capture much of the physics of 
multi-species membrane junction separation induced by adhesion. I show that not only the 
difference of junction height, but also the difference of junction rigidities, and the membrane- 
mediated interactions between the junctions play important roles in the junction separation. 
The fluctuation analysis also shows that current analysis does not provide the complete phys- 
ical picture for systems with very soft or very short junctions, or in the situation when the 
junction densities are extremely low, where the thermally activated membrane fluctuations 
or the Helfrich repulsion between the membranes become important interactions in the sys- 
tem [I9|. In this regime, numerical simulations |7], |2D| should provide valuable information 
on the distribution of the junctions, as well as a complete picture of the phase diagram which 
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includes the bind/unbinding transition between the membranes, and adhesion- induced junc- 
tion separations. 
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Appendix 

In this Appendix, I discuss the details of some calculations mentioned in the text. For 
simplicity I define 

5m + (r) = Ai50i(r) + A 2 50 2 (r) 

5m_(r) = ArfMr) - A 2 6<h(r). (31) 

First, an integral which is very useful for the rest of this Appendix is calculated. 

d 2 q A(q)B(-q) f 2 f 2 f d 2 q A{v)B{v')e^-^ 
(2tt) 2 q 4 + m + J J T J (2vr) 2 g 4 + m+ 

f ,o f , 9 , 1 f°° dx xJ (x\r — r'lmV 4 ) „. , „. 

= I dr r r ^ml a? x* + i - 4(r)B(r) 

= f d 2 r [ d 2 r'G{\T-v'\m 1 l A )A(v)B{r'), (32) 



where G(x) is a MeijerG function fll8"| . Also. 



7T, X <C 1 

0, x > 5. 



The shape of G(x) is plotted in Fig. 6. 

Now I consider the nonlocal interaction between the junctions in the Gaussian approxi- 
mation. Neglecting the first term of Eq. fl22|) , the second term can be expressed by 

\l M Sm + (q) + A fe (5m_(q)| 2 
G ~ V g 4 + rg 2 + m+ 

ci 2 g |/ M <5m + (q) + A fe £m„(q)| 2 
(2vr) 2 g 4 + Tq 2 + m+ 
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- / d 2 r / d 2 r' 



^LG(|r-r'|mf) 



n? 



Ai^i(r) 



+ , 



1 + - A*(r) 
m+7 



1- — ) Ai^i(r') - 1 + — A 2 ^ 2 (r') 

m + / V 777,4- / 



(33) 



where the last expression holds when T < y/m + , and the integral in Eq. ( p2|) is used to 
calculate the Fourier transformation from the momentum space to the real space. The 
range of this membrane mediated interaction between the junctions is determined by the 
shape of G(x), which sets the length scale of this interaction to m+ . 

Next I show that the interaction between the junctions due to the contribution of Hi is 
of the form in Eq. (p5|). The effective interaction free energy between the junctions is 



-In 



D[h]e 



-Hm-Ho-H^-H. 



H m + H 6 -ln 



D[5l]e 



-H -H! 



(34) 



When the contribution of Hi is included by a one-loop calculation, 

F c = H M + H (/) + F G + < Hi > -- (< Hi > -<Hi > 2 
= Fm + Fa + Fi oop + H<p. 



Here 



F M = H M , F G = - ln[ / D[5l}e- H °], and < O > ~- 



JD[5l] Oe~ H ° 



(35) 



(36) 



JD[6l\e- B ° 

for any O. The calculation of < Hi > is straightforward, which in the case Y < v /m + leads 
to 

< Hi > = l - J d 2 r{KMi{v) + A 2 5</> 2 (r)) < (5l(r)) 2 > 



d 2 q 



(2tt) 2 q A + iy + m+ 



(A 1 50 1 (r)+A 2 50 2 (r)) 



16, /mZ 



dMAi&MO + A^r)). 



(37) 



Terms of higher order than 5(j) 2 a have been neglected. The calculation for < Hi > is longer 
but also straightforward, 

<Hl> = - J d 2 r J d 2 r' < 8l(r)5l(r)5l(r')5l(r') > 5m + {r)5m + {r') 

1 f d 2 q r d 2 q' 1 1 . . . , l2 

' ' |5m + (q)| 2 + < Hi 



2 J (2tt) 2 J (27r)V 4 + iY 2 + m + (q + q') 4 + r(q + q') 2 + m + 
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Again, terms of higher order than b~<\) 2 a are neglected. When the contribution from the surface 
tension can be neglected, the following form provides a good approximation of < H\ >o — < 
Hi >q, |2T] this form gives the correct asymptotic behavior of the true result at large and 
small q limits, 

<H(> Q -<Hi> 2 « — — / dy 



i= / d 2 r f d 2 /G?(|r-r / |(4m+) 1 / 4 ) (Ai<tyi(r) + A 2 ^ 2 (r)) 



x(A 1 ^ 1 (r / ) + A 2 50 2 (r / )) (39) 

Putting < iJi > and < H\ >o — < H\ >\ together leads to the resulting expression of 
F c in the one- loop order, which is given in Eq. ( p4|) and Eq. (|25|). The real space form of 
< H\ >o — < Hi >q also shows that the membrane fluctuation induced interaction between 
the junctions is short ranged with a characteristic length on the order of m^ 1 ^ 4 
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Figure Captions 
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• Fig.l Schematics of the system. The membrane heights are Zi(r) and z<i (r) from the 
reference plane. There are two types of receptors in one membrane, two types of 
ligands in another membrane. Two types of junctions can be formed from the ligands 
and receptors. They have different natural lengths hi and h 2 . In general, different 
types of junctions also have different rigidities. The softer junctions are easier to be 
stretched or compressed from their natural length. 
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• Fig.2 The shape of g(<f>) with different values of //_ when A = 0.2, and A|A + = 1.998. 
Solid line: \x = —0.404, dashed line: \x = —0.4045, dash-dotted line: /x = —0.405. 
Phase coexistence occurs at /i ~ —0.4045 even though A^A + < 2.0. 
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Fig.3 g(4>) in the phase coexistence for different values of A with A^A + = 2.04. Solid 
line: A = 0, dashed line: A = 0.025, dash-dotted line: A = 0.05. A = curve is 
symmetric around <fi = 0, A > curves shows that the positions of the minima are 
shifted towards smaller values, i.e., softer junctions are the favored species. 
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• Fig.4 Schematics of the phase coexistence curves near A|A + = 2 for A <^ 1 for A = 
(thick solid line), A = ±Ai (thick short-dashed lines), and A = ±A 2 (thick long-dashed 
lines). Ai > A 2 > 0. The thin dashed curve is the position of the end points of the 
phase coexistence curves, this is given by A|A + 2(1 — 9A 2 /4), /i_ « — 2A. The 
curves move towards the A = phase boundary as A|A + increases because the effect 
of junction height mismatch becomes more important. 
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• Fig. 5 Membrane- mediated interaction revealed by the Gaussian approximation. This 
interaction comes from the bilinear coupling between SI and S(f> a . (A) A small region 
which has higher density in the junctions with greater natural height (or lower density 
in the junctions with smaller natural height) induces a positive 81. Two regions with 
positive SI can reduce the bending elastic energy of the membranes by moving close to 
each other. Similarly, a region with negative SI attracts another region with negative 
SI due to the cost of membrane bending energy. (B) A small region with positive SI 
repels with a region with negative SI because of the bending elastic energy cost of the 
the high curvature region between them. 

G(x) 
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• Fig.6 The shape of the MeijerG function G(x). Although G(x) oscillates very weakly, 
and has a local minimum close to x — 5, the important feature of G(x) is that this 
function is vanishingly small when x > 5. 
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